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We develop a new ab initio many-body approach capable of describing simultaneously both bound 
and scattering states in light nuclei, by combining the resonating-group method with the use of realis- 
tic interactions, and a microscopic and consistent description of the nucleon clusters. This approach 
preserves translational symmetry and Pauli principle. We present phase shifts for neutron scatter- 
ing on ^H, "^He and ^"Be and proton scattering on '^'^He, using realistic nucleon-nucleon potentials. 
Our A — 4 scattering results are compared to earlier ab initio calculations. We demonstrate that a 
proper treatment of the coupling to the n-^^Be continuum is essential to explain the parity-inverted 
ground state in ^^Be. 
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The development of a quantitative microscopic the- 
ory of low-energy nuclear reactions on light nuclei is 
key to further refine our understanding of the funda- 
mental nuclear interactions among the constituent nu- 
cleons, providing, at the same time, accurate predictions 
of crucial reaction rates for nuclear astrophysics. The 
enormous difficulties involved with such a project are 
immediately apparent if one considers the nearly-total 
lack of ab initio calculations for scatteringprocesses in- 
volving more than four nucleons overall Neverthe- 
less, a breakthrough is within reach today thanks to 
the significant progress achieved in the ab initio descrip- 
tion of the structure of light nuclei and the advent of 
modern high-performance computers. In this Letter we 
combine a microscopic-cluster technique, the resonating- 
group method (RGM) 0, and a very successful structure 
approach, the ab initio no-core shell model (NCSM) 3], 
into a new microscopic theory {ab initio NCSM/RGM) 
capable of treating bound and scattering states of light 
nuclei in a unified formalism, starting from the funda- 
mental inter-nucleon interactions. Within this new ap- 
proach we study the n-'^H, n-^He, n-^°Be, and p-^'^^He 
scattering processes, and address the parity inversion 
of the ^^Be ground state (g.s.), using realistic nucleon- 
nucleon (NN) potentials. 

We start from the wave function for a scattering pro- 
cess involving pairs of nuclei that can be cast in the form 



(1) 



through an expansion over binary-cluster channel-states 
of total angular momentum J, parity tt, and isospin T, 
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The wave functions of the {A — a)- and a-nucleon clus- 
ters are antisymmetric under exchange of internal nu- 
cleons, and depend on translationally-invariant internal 



coordinates. They are eigenstates of the 



(A-a) 



and 



intrinsic Hamiltonians with spin, parity, isopsin and 
additional quantum numbers Ii,T:i,Ti, and a^, respec- 
tively, where i — 1,2. The clusters centers of mass are 
separated by the relative vector r^-a.a- Relative angu- 
lar momentum and channel spin are denoted by i and 
s, respectively. The inter-cluster anti-symmetrizer for 
the {A — a, a) partition in Eq. ([1]) can be schematically 
written as A = [{A - a)!a!/^!]i/2 ^p(-)PP, where P 
are permutations among nucleons pertaining to different 
clusters, and p the number of interchanges characterizing 
them. The coefficients of the expansion with respect to 
the channel index ^ — {A~aaiIi^Ti] aa2l2^T2] sC\ are 
the relative-motion wave functions (r) , which repre- 
sent the unknowns of the problem. They can be deter- 
mined by solving the many-body Schrodinger equation in 
the Hilbert space spanned by the basis states Ai, l^lr'^), 



where the integration kernel is given by: 
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Here E is the total energy in the center-of-mass (cm.) 
frame, and H is the intrinsic j4-nucleon microscopic 
Hamiltonian, which it is useful to decompose into, e.g.: 



H 
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Further, T,.g\{r) is the relative kinetic energy and Vrci 
is the sum of all interactions between nucleons belong- 
ing to different clusters after subtraction of the avarage 
Coulomb interaction between them, explicitly singled out 
in the term Vc{r) = Zi,yZ2i^e^ /r, where Zi^, and Z21, are 
the charge numbers of the clusters in channel ly. 

We obtain the cluster eigenstates entering Eq. ([2]) 
by diagonalizing Hi^A-a) and iJ(a) in the model space 
spanned by the NCSM basis. This is a complete har- 
monic oscillator (HO) basis, the size of which is defined 
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TABLE I: Calculated and ^He g.s. energies (in MeV), n-'^H, n-*He and p-^He phase shifts (in degrees), and n-^H and 
n-^He total cross sections (in barns) for increasing A'max at HQ, = 18 MeV, obtained using the Viowk NN potential (derived 
from AV18 with cutoff A = 2.1 fm~^) [||. Only the g.s. of the '^H and ^He nuclei were included in the scattering calculations. 



by the maximum number, A^max, of HO quanta above 
the lowest configuration shared by the nucleons. Thanks 
to the unique properties of the HO basis, we can make 
use of Jacobi-coordinate wave functions [J] for both nu- 
clei or only for the lightest of the pair (typically a < 4), 
and still preserve the translational invariance of the prob- 
lem. In the second case we expand the heavier cluster on 
a Slater-determinant (SD) basis, and remove completely 
the spurious cm. components in a similar fashion as in 
Refs. [^, @] . We exploited this dual approach to verify our 
results. The use of the SD basis is computationally ad- 
vantageous and allows us to explore reactions involving p- 
shell nuclei. In calculating ([4]), all "direct" terms arising 
from the identical permutations in both Au and Ay' are 
treated exactly with the exception of ($;^vT| Vrei 
The latter and all remaining terms are obtained by ex- 
panding the Dirac 5 of Eq. ([2]) on a set of HO radial 
wave functions with identical frequency f2, and model- 
space size iVinax consistent with those used for the two 
clusters. In this respect we note that Vroi is localized also 
in presence of the Coulomb force. We solve Eq. Q by 
means of the coupled-channel i?-niatrix method on a La- 
grange mesh f?] imposing either bound-state or scattering 
boundary conditions for -^(r) at large r. 

All calculations in the present paper were carried out 
using binary-cluster channels ([2|) with a = 1. We first dis- 
cuss results obtained limiting the expansion ([T|) to config- 
urations with the {A — l)-cluster in its g.s. Table |T] shows 
the behavior with respect to TVmax of selected A = 4, 5 
data obtained using the Viowk NN potential Q . A satis- 
factory convergence of both g.s. energies and scattering 
data is reached starting from A^max = 17, for the four-, 
and A'niax = 15 for the five-nucleon systems, respectively. 



In what follows we present results obtained using ef- 
fective interactions derived from the underlying realis- 
tic NN potential. Vat, through a unitary transformation. 
Starting from the relevant two-nucleon Hamiltonian (for 
notation and definitions see Ref. — H02 + V12, 

with V12 — VN{\/2r) — 'mQ'^r^/A, the cluster eigenstates 
are obtained employing the usual NCSM two-body eflFec- 
tive interaction V2cS — H2cS — Hq2, where i?2cff is the 
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FIG. 1: (Color online.) Calculated phase shifts for n-^H scat- 
tering as a function of the relative kinetic energy in the cm. 
frame Skin, using the N^LO NN potential [§]. Only the g.s. 
of the {A — l)-cluster was included in the present calculation. 
Dependence on the model-space truncation A'max at hO. — 22 
MeV compared to AGS resuhs of Ref. [13, (H- 
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FIG. 2: (Color online.) Same as Fig. [l]for p -^He scattering. 



Hermitean effective Hamiltonian. However, in place of 
the bare NN potential entering Vid we adopted the new 
effective interaction Vjgg = ^?2eff ~ ^2cff ' where H!^^^ is 
the effective Hamiltonian derived from H2 ' = ^02 + ^12 > 
with = -mOF-r'^/A. Note that V^ff V^r in 
the limit A'^^ax 00. Figures [H [2] and the right 
panel of Fig. [3] present A = A,5 scattering phase shifts 
for the high-quality NN potential derived within chiral 
effective-field theory at next-to-next-to-next-to-leading 
order (N^LO) [9]. For the whole energy range, we find 
less than 2 deg absolute difference between the phases 
obtained in the largest and next-to-largest model spaces, 
a sign of convergence. The only exception is represented 
by the ^^3/2 phase shifts of the n-a system, for which this 
difference rises up to 5 deg in the range 1 MeV< i?kin < 4 
MeV. As a comparison, we show in the left panel of Fig. [3] 
the n-a phase shifts obtained with the (bare) Viowk in- 
teraction. The convergence rate is clearly much faster. 

In order to verify our approach, in Fig. [1] and [2] we 
compare our n-^H and p-^He results to earlier ab ini- 
tio calculations performed in the framework of the Alt, 
Grassberger and Sandhas (AGS) equations 13, U], using 
the same N'^LO NN potential. We note that in general 
the agreement between the two calculations worsens as 
the relative kinetic energy in the cm. frame, -Ekin, in- 
creases. For the P-waves in particular we can reasonably 
reproduce the AGS calculation for energies within 1 MeV 
off threshold, while we can find differences as large as 17 
deg ('^P2) at Ekin — 2.6 MeV. These discrepancies are 
due to the influence, increasing with energy, played by 
closed channels not included in our calculations, such as 
those with the A — 1 = 3 eigenstates above the I^'^ =5^ 
g.s., and (A — a = 2, a — 2) configurations, present in the 
AGS results. In Ref. [13] it was shown that the omission 
of three-nucleon partial waves with ^ < /i < | leads to 
effects of comparable magnitude on the AGS results at 
£;kin = 3 MeV, especially for the ^Si,^ Pi and ^P2. 

We explore the effect of the inclusion of higher excited 
states of the {A — l)-cluster on the n-a scattering phase 
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FIG. 3: (Color online.) Same as Fig. [T] for n-a scattering 
with Viowk NN potential [1] at fifJ = 18 MeV (left panel), and 
N^LO NN potential ^ at HQ = 19 MeV (right panel). 



shifts. Channels with a > 1 have here a much suppressed 
effect due to the large binding energy of the ^He nucleus. 
Figure 2] shows the influence of the six lowest excited 
states of "'He on the n-a phase shifts. The I^^Ti = 0+0 
excited state affects only minimally the ^S'1/2, leaving the 
P phase shifts unaltered. On the contrary, we find larger 
deviations on the ^Pi/2 and ^^3/2 phase shifts, after the 
inclusion of the 0^0, 1~0, and 1~1 states for the first, 
and of the 2~0 and 2~1 states for the second. These 
negative-parity states do not influence the ^51/2. 

In Fig. [5] the n- and p-a phase shifts obtained with 
the N'^LO NN potential, including the first six '^He ex- 
cited states, are compared to the results of an accurate 
multi-channel i?-matrix analysis of the nucleon-a scatter- 
ing data [13]. The ^S'1/2 phase shifts are in good agree- 
ment with experiment, also in presence of the Coulomb 
repulsion between proton and a particle. The magnitude 
of the ^1)3/2 phase shifts is also qualitatively reproduced. 
On the contrary, the P phase shifts present both insuffi- 
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FIG. 4: (Color online.) Influence of the lowest six exited 
states (0+0, 0"0, 1"0, 1" 1, 2"0, 2" 1) of the a particle on the 
n-a phase-shift results for the N^LO NN potential j9|. 
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FIG. 6: (Color online.) Calculated 1/2 n-^^Be phase shifts 
as a function of Skin, using the CD-Bonn NN potential [l^ 
at hO. = 13 MeV. NCSM/RGM calculation as in Tab.|lll 



FIG. 5: (Color online.) Calculated n-a (left panel) and 
a (right panel) phase shifts for the N^LO NN potential 
including the "He g.s., 0+0, 0^0, 1^0, 1^1, 2~0, and 2" 
states, compared to an _R- matrix analysis of data (+) [l^. 



cient magnitude and splitting with respect to the predic- 

1 + 

tions of the i?-niatrix analysis. The ^ channel is dom- 
inated by the repulsion between nucleon and a particle 
induced by the Pauli exclusion principle. Consequently, 
the short-range details of the nuclear interaction play a 
minor role on the ^>5'i/2 phase shifts, for which, as shown 
in Fig. [21 wc find very similar results using the Viowk po- 
tential. On the other hand, the latter figure shows also 
that the and ^^3/2 phase shifts are sensitive to the 

interaction model, and in particular to the strength of 
the spin-orbit force. The present discrepancy with re- 
spect to experiment is due to the omission of the three- 
nucleon terms of the chiral interaction at order N^LO, 
which would lead to an enhanced spin-orbit splitting. 

To show the promise and flexibility of our approach, we 
present in Tab.Hlland Fig.[6]results for a much heavier {A 
= 11) system. The parity-inverted g.s. of ^^Be, one of the 
best examples of disappearance of the N — 8 magic num- 
ber with increasing N/Z ratio, was so far left unexplained 
by ab initio calculations [3| • The HO asymptotic behav- 
ior of the ^^Be wave function in the standard NCSM does 







lOBe 


"Be(i-) 


"Be(i+) 






Eg.s. 


E 


Eth 


E 


Eth 


NCSM [13, JJ 


8/9 


-57.06 


-56.95 


0.11 


-54.26 


2.80 


NCSM [13,24]," 


6/7 


-57.17 


-57.51 


-0.34 


-54.39 


2.78 


NCSM/RGM" 






-57.59 


-0.42 


-57.85 


-0.68 


Expt. 




-64.98 


-65.16 


-0.18 


-65.48 


-0.50 



"present calculation 

TABLE II: Calculated energies (in MeV) of the ^"Be g.s. and 
of the lowest negative- and positive-parity states in ^^Be, ob- 
tained using the CD-Bonn NN potential [l^ at ftfi = 13 MeV. 
The NCSM/RGM results were obtained using n-|-^°Be config- 
urations with A^max = 6 g.s., 2+, 22 , and 1+ states of ^"Be. 



not favor extended n-^'^Be configurations, thus enhanc- 
ing the relative kinetic energy repulsion, and preventing 
the experimentally-observed inversion between i and 



states. Using the CD-Bonn NN potential 15|, we ob- 



serve a dramatic (~ 3.5 MeV) increase in the Be 



1 + 
2 

state binding energy leading to the g.s. parity inversion, 
when the n-^'^Be relative motion is treated within the ab 
initio NCSM/RGM approach. 
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